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1. Introduction 

HERWIG^ g is a general-purpose event generator for high-energy processes. A major 
new feature is the inclusion of supersymmetric (SUSY) processes, made available for the 
first time in HERWIG 6.1 and extended in subsequent versions. This article describes 
the implementation, defines the parameter and mixing conventions, and lists the hard 

^The latest version of the program, together with all relevant information, is always available via the 
official HERWIG page: ^ittp : / /hepwww . rl . ac .uk/theory/ seymour/herwig/|. 
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scattering matrix elements. For the details of the input file format and for a more general 
description of our implementation, as well as the description of the Standard Model (SM) 
physics implemented in the generator, the reader should consult Ref. |Q. This article 
supplements that paper and extends it to version 6.4. The exact version number of the 
code described in this article is 6.400. 

The key features of the HERWIG 6.4 supersymmetry implementation are as follows: 

• A general Minimal Super symmetric Standard Model (MSSM) particle spectrum. 

• Both hadron-hadron and lepton-lepton collisions. 

• MSSM 2 — > 2 sparticle production subprocesses. 

• 2^1, 2^2 and 2^3 Two Higgs Doublet Model (2HDM) Higgs boson produc- 
tion subprocesses in hadron-hadron collisions, and 2 — > 2 Higgs boson production 
processes in lepton-lepton collisions. 

• All R-parity violating 2 — > 2 scattering subprocesses in hadron-hadron collisions 
that can, but do not necessarily, have an s-channel resonance, and 2^2 scattering 
processes in lepton-lepton collisions, not necessarily with an s-channel resonance. 

• Cascade decay of heavy objects according to branching fractions taken from an input 
data file, incorporating 3-body and 4-body matrix elements. 

• Colour coherence between emitted partons. 

• Spin correlations in 2 ^ 2 R-parity conserving production processes and all cascade 
decay processes. 

• Beam polarisation in lepton-lepton collisions. 

The masses, total decay widths and branching fractions, the mixings of gauginos, 
third generation sfermions and Higgs bosons, the bilinear, trilinear and R-parity violating 
trilinear couplings, are read in from a data file. A separate program (ISAWIG) is available 
(from the HERWIG web page) to convert the data from ISAJET g and HDECAY g to a 
form suitable for HERWIG and modify it to include the R-parity violating sector. 

In our implementation, super symmetric particles do not radiate gluons. This assump- 
tion is reasonable if the decay lifetimes of the coloured sparticles (spartons) are much 
shorter than the QCD confinement scale. In particular, a stable gluino Lightest Supersym- 
metric Particle (LSP) is not simulated. The simulation of a long-lived light top squark 
would neglect effects due to gluon emission and hadronization. 

CP-violating phases are not included for either the mixings or the couplings. 

R-parity violating lepton-gaugino and slepton-Higgs mixing is not considered. 
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2. Parameter and mixing conventions 

There are several different conventions that are commonly adopted in the literature for 
the MSSM Lagrangian. While most authors have chosen to follow the conventions used 
in Refs. |5|, other conventions, for example Ref. [Q, are also commonly used. The 
conventions we have used in general follow those of Refs. |Q, ^ • In this section we present our 
conventions for the mixing matrices for the charginos, neutralinos and sfermions, together 
with the conversion between the conventions of Refs. ^ and Ref. [^. 

2.1 Particle content 

The MSSM particle content, expressed in terms of its mass eigenstates, is given in Tab. 1. 



spin 


spin 1/2 


spin 1 


down-type squarks 
dfl, Si, Sfl, 61,62) 


down-type quarks 
{d,s,b) 




up-type squarks 


up-type quarks 
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charged sleptons 
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(e-,^-,r-) 
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(l/e,t'/,,i^T) 
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9 
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9 


neutral Higgs 
(/^o,i70,v4") 


neutralinos 

X° = 1 - 4) 


photon and 


charged Higgs 
H+ 


charginos 
X+ = 1 - 2) 


W+ 




Goldstino 
G 





Table 1: Particle content of the MSSM, expressed in terms of its mass eigenstates. The 
Goldstino G is the spin-1/2 component of the gravitino. We neglect mixing for the first 
two generation sfermions. 

In defining our conventions we need to write the particle content in terms of the 
interaction eigenstates. This is shown in Tab. 2. Superscripts c refer to charge conjugation, 
which, stated in the four-component notation, is defined by: 

= Ci^T ^ ^^2^0^T_ ^2.1) 

We adopt the chirality basis. The two-component case, which we refer to when dis- 
cussing mixing conventions, follows trivially. Acting on chirality-left and chirality-right 
two-component spinors, we have: 

{^lT = i^^^l = ^ij^jL^ (2-2) 
{cPnr = -i^^(P*R = -eiA- (2-3) 
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The Levi-Civita tensor eij is defined by £12 = —£21 = l^eii = £22 = 0. 
2.2 Parameter conventions 

In terms of the scalar fields defined in Tab. 2, we write our scalar superpotential as follows: 



(2.4) 



The above definition is consistent with Refs. (Bj, The sign of the term proportional to the 
up-type quark coupling Yukawa matrix fu is such that the Yukawa couplings are positive. 
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Table 2: Field content of the MSSM, expressed in terms of its interaction eigenstates. 
Generation and colour indices are suppressed. 



For the R-parity violating scalar superpotential, we take: 



f^aL^H^ + -XabcLaLlEc + X'abcLaQlDc 



+ 2'^abc^ciC2C3Ua^ D'^^ ■ (2.5) 



a,b,c are the generation indices. In the last term, 01,02,03 are the triplet colour indices 
and the totally antisymmetric tensor is defined by £123 = 1- Summation over colour indices 
is implicit in the other terms. We omit the bilinear term proportional to k^. 
We introduce our soft SUSY-breaking bilinear and trilinear terms as follows: 



Vi 



fiBHlH^^ + feArHiVE + foAbHlQW - fuAtW^Q^U + (h.c 



(2.6) 



This is consistent with the definition in Ref. We neglect A term contributions to the 
first two generations. The contribution of the potential to the Lagrangian is £ = —V in 
general, (h.c.) stands for Hermitian conjugation. 
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Our soft SUSY-breaking masses are defined schematically as follows: 

Vm = Mlct^H + M^^^, (2.7) 

where cp are the sfermion fields and ip are the gaugino and gluino fields in the fom'- 
component notation. Summations over colour and weak isospin indices are implicit. We 
write the gaugino soft-breaking masses as = Mi and = M2. In the notation of 
Refs. IQ, ^, M' and M are used instead of Mi and M2, respectively. However, the notation 
Ml and M2 has become more common. 

The above conventions for the signs of the and the A terms are the same as those 
adopted in ISAJET and Haber and Kane This also agrees with the conventions 
of the other major SUSY event generators, SPYTHIA §] and SUSYGEN [§. While these 
conventions are the same internally, ISAJET uses a difi^erent notation which corresponds to 
the following: 

2mi = -ij, (2.8) 
m = -Ml, (2.9) 
^12 = -M2, (2.10) 

where mi, /ii and ^2 are in the notation of ISAJET. 

2.3 Mixing conventions 

While internally HERWIG uses the conventions of Refs. |Q, ^, Q, the mass and decay spectra 
are fed in from a data file. In particular, we have provided a separate code (ISAWIG) for the 
conversion of data from ISAJET |^ . This way, the masses of the sparticles and their R-parity 
conserving decay rates are calculated using ISAJET and the R-parity violating section is 
added by ISAWIG. ISAJET also calculates the mixing matrices for the electroweak gauginos 
and scalar fermions. ISAJET uses the conventions of Refs. (2|, |l^ and it is essential that 
the conversion between the two formalisms is done correctly. 

We first consider the conversion between the two formalisms for the gauginos and then 
for the left -right sfermion mixing. We compare the Lagrangian of Ref. Q against that of 
Refs. ^ to obtain the conversions between them. We adopt the conventions of Refs. |5| 
for the chargino and neutralino mixing matrices. 

2.3.1 Charginos 

In the notation of Ref. ||^], we define the following doublets of two-component spinors for 
the charginos in their interaction eigenstates: 

^P+ = (2.11) 

= [-i\-,i'H,)^ (2.12) 

where are the charged Winos, tlj^_^ is the charged Higgsino associated with the Higgs 
that gives mass to the down-type quarks and i/^^jj is the charged Higgsino associated with 
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the Higgs doublet that gives mass to the up-type quarks. The mass term in the Lagrangian 
from eqn. (A. 2) of Ref. is: 

^chargino = -\ (V+ V'") ^ ^ ^ ^ ^- ^ + (^-C-), (2-13) 

where the mass matrix is given by: 

X={ MH.V2sin/3Y 

In Ref. 1^] this is diagonahzed in the two-component notation by defining the mass eigen- 
states: 

Xt = V^,^+, (2.15) 
xr = U^J^J, (2-16) 

where U and V are unitary matrices chosen such that 

U*XV-^ = Md. (2.17) 

Md is the diagonal chargino mass matrix. The four-component mass eigenstates, the 
charginos, are defined in terms of the two-component fields as: 

We need to express the chargino mass terms in the Lagrangian in a four-component notation 
so that we can compare this Lagrangian with the conventions of Ref. |^. We define the 
four-component spinors as in Ref. [^: 

Using this notation, we can express the Lagrangian in the four-component notation: 



>Cchargino = " {w+ H+) [XPl + Pr] ^+ j ' (2-20) 

The chirality projection operators are defined as Pl/r = (l^/+75)/2- Eqn. ( |2.20| ) has 
a similar form to the chargino mass term^ in eqn. (2.1) of Ref. Q: 

-Cchargino = " (A x) [M.harge^'L + Afi,,g,Pfi] ( M , (2.21) 



^This Lagrangian is taken from Ref. which corrects a sign error in the off-diagonal terms in Ref. 
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where 



charge 



/i2 -gv 
—gv 2mi 



M2 V2Mw cos P 
\/2Mw sin P /i 



(2.22) 



where v and v' are the vacuum expectation values for the Higgs fields that give mass to 
the up- and down-type quarks, respectively. We now come to the problem of comparing 
the notation of Ref. against that of Ref. Q. The Wino and Higgsino fields of Ref. Q 
are the charge conjugates of those used in Ref. |||], so that the following transformation is 
required: 



This gives the chargino mass matrix in the form: 

-^charge ~^ charge 

where 



•^char 



gino 



W+,H+ 




charge 





gv 


=( 


I gv' 


-2mi . 





M2 V2Mwsm(3 
V2Mw cos /3 /X 



(2.23) 
(2.24) 



(2.25) 



(2.26) 



This agrees with eqn. (2.20) apart from the overall sign. We need to express the fields in 
terms of the ISAJET mixing matrices. The ISAJET mixing matrices are given in eqns. (2.10) 
and (2.11) of Ref. @ as: 




cos 7L 



;Sin7L \ / A 



i-iy+w+ 



sm7L COS7L 



9y cos 'jR -Oysin-^R \ ( 
sin jR COS jR 




(2.27) 
(2.28) 



The mixing angles 7^ and jr, and the sign functions 6x, Oy, 6+, and 9- are defined in 
Ref. i. 

We can transform these equations into the notation of Ref. 0] with the identification: 



This gives: 




xt = wi, 
xt = wi. 



-sin7i? -OyCosjR 
cosjR OySin-fR 



cos 7l 



-f^x COS7L 

9^ sin 7L 




(2.29) 
(2.30) 



(2.31) 
(2.32) 
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These mixing matrices are now in a form that allows us to compare them with the notation 
of Ref. §, eqn. (A.13): 




Pr 



Uu U21 
U12 U22 




(2.33) 
(2.34) 



By comparing eqns. ( p. 32 ) and ( 2.34| ) we obtain the mixing matrices in the notation of 
Ref. 0: 



U 



V 



- sm jL 

COS 7L 



— COS 7l 
9x sin7L 



-sin7K(-l)''- -cos7r(-1)^+ 
-0j,cos7r(-1)^- eySm-fR{-lf+ 



(2.35) 
(2.36) 



It should be noted that we adopt the opposite sign convention for the chargino masses due 
to the sign differences in the two Lagrangians. 

2.3.2 Neutralinos 

We define a quadruplet of two-component fermion fields for the neutralinos in the interac- 
tion eigenstate as: 

V'° = HA',-a3,V'?,,,V'Hj, (2.37) 

where A' is the Bino, is the neutral Wino, ip^^ is the Higgsino corresponding to the 
neutral component of the Higgs doublet giving mass to the down-type quarks and ipjj^ is 
the Higgsino corresponding to the neutral component of the Higgs doublet giving mass to 
the up- type quarks. The Lagrangian for the neutralino masses from equation (A. 18) of 
Ref. H is: 



C 



1 



neutralino 



i(/)'"yv° + (h.c.), 



(2.38) 



where 



Y 



Ml 


-Mz sin 6]y cos /? 





M2 

Mz cos 9w cos P 



\ Mz sin 6w sin /3 —Mz cos 6w sin /3 



-Mz sin 9w cos /3 
Mz cos 6w cos P 




Mz sin 6w sin /? \ 
-Mz cos 6w sin (3 





(2.39) 

In Ref. 1^ , the Lagrangian is diagonalized in this two-component notation to give the mass 
eigenstates. The diagonahzation is performed by defining two-component fields: 



(2.40) 



where is a unitary matrix satisfying N* Y N ^ = Nd , and Nd is the diagonal neutralino 
mass matrix. The four-component mass eigenstates, the neutralinos, can be defined in 
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terms of the two-component fields: 




(2.41) 



Rather than adopting this approach we need to express the neutrahno mass terms in the 
four-component notation before performing the diagonahzation in order to compare this 
Lagrangian with that of Ref. |Q| . We use the standard procedure of Ref . [|[ to express this 
in the four-component notation by defining four-component Majorana fields: 



B 





Hi 



Ho 



H2 



(2.42) 



H2^ 



As defined in Tab. 2, B is the Bino field, W3 is the neutral Wino field. Hi is the field for 
the Higgsino associated with the neutral component of the Higgs doublet that gives mass 
to the down-type quarks and H2 is the field for the Higgsino associated with the neutral 
component of the Higgs doublet that gives mass to the up-type quarks. 
This gives the Lagrangian in the four-component notation: 



/:neutralino = (i? ^^3 ^l,i^2) [Y Pl + Y Pr] 



( B\ 

Hi 
\H2j 



(2.43) 



We can now compare this with the relevant Lagrangian, in the four-component notation, 
given in Ref. H. This Lagrangian^ is from Ref. S, eqn. (2.2): 



neutralino 



■, h °, A3, Ao ) [Mneutral-Pi + neutral -Pr] 



h 



A3 
VAo/ 



(2.44) 



/i° is the Hig gsino partner of the neutral component of the Higgs doublet that gives mass 
to the up- type quarks, h'^ is the Higgsino partner of the Higgs boson that gives mass to 
the down-type quarks, A3 is the neutral Wino and Aq is the Bino. The mass matrix is given 
by: 



-^neutral 





—2mi 
-gv/V2 
g'v/V2 



—2mi 


gv'/V2 
-g'v'/V2 



-gv/V2 

gv'/V2 





g'v/V2 

-g'v'/V2 





\ 



(2.45) 



^This Lagrangian is taken from Ref. U/w which corrects a sign error in Ref. 
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It is easier to compare this with eqn. ( ^.43 ) after reordering the entries and reexpressing 
it in terms of Mz-, (3 and 9w- This gives: 



/Ao\ 



-Cneutralino — - ( Aq, A3, /i °, /l° ] [M^^^^^^iPl + M^^^^^j.^iPji] 



(2.46) 



where: 



M' 

neutral 



/ Ml 


Mz sin 9w cos /? 
\—Mz sin 9\Y sin /3 





M2 

-Mz cos 9\Y cos /3 
cos 9w sin /? 



Mz sin 9w cos /3 
-M^ cos 9w cos /3 




-M^ sin sin /? \ 

M^ cos sin P 

-/^ 




(2.47) 

This equation then agrees with eqn. ( |2.43| ) up to an overall sign, provided we make the 
following identification: 



B 

Hi 
H2 



A3, 



(2.48) 
(2.49) 
(2.50) 
(2.51) 



The convention from ISAJET for the mixing, taken from eqn. (2.12) of Ref. [0], is 



^^3 



Z3 
Z4J 



/ (1) (1) (1) (1) \ 

v) vX v-, V A \ 



(2) ^,(2) ^,(2) ^,(2) 



Vn 



Vn 



Va 



(3) „,{3) .,(3) .,(3) 
1 

(4) 
1 



,(4) „(4) ,,(4) 



A3 

VAo/ 



(2.52) 



where Zj are the mass eigenstates obtained by diagonalizing the mass matrix given in 
eqn. ( 2.45| ), uj*^ are the elements of the mixing matrix, and 9i is zero (one) if the mass of 
Zi is positive (negative). 

After reordering, we get: 



(-^75)'^ Z2 
i-il^f' Zj 

\{-ii^tz^) 



V 



(2) ^,(2) ^,(2) ^,(2) 



V 



V 



(3) .,(3) .,(3) „,{3) 



V 



(4) ,,(4) ^,(4) ^ (4) 



W 
Hi 
\H2j 



(2.53) 



We can therefore obtain the mixing matrix in the notation of Ref. by making the 
identification: 



iV^l 



^2 > 



(2.54) 
(2.55) 
(2.56) 
(2.57) 
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Again, we need to adopt the opposite sign convention for the neutrahno masses. 



2.3.3 Left right sfermion mixing 

In addition to the mixing of the neutraUnos and charginos, we need to consider the left- 
right mixing of the sfermions. In general, as the off-diagonal terms in the mass matrices are 
proportional to the fermion mass, these effects are only important for the third generation 
sfermions, i.e. stop, sbottom and stau. We describe our convention for the top squarks. 
The others are treated analogously. 

The fohowing mass matrix for the top squarks uses the conventions^ of Ref. |^ and is 
taken from eqn. (4.17) therein. 

/^^| + ^|cos2/?(i-|sin2 0H/)+m2 mt (A - cot P) \ 

* \^ mt{At- fi cot P) Mi + |m| cos 2/3 sin^ 9w + j 

where Mg and are soft SUSY-breaking masses for the left and right top squarks, 
respectively. At is the trilinear soft SUSY-breaking term for the interaction of the left and 
right stop squarks with the Higgs boson. This compares with the ISAJET matrix from 
Ref. pi: 



+ M| cos 2/3 (i - I sin^ Ow) + mj -rrit {At - fi cot /3) 

+ |m| cos 2/3 sin^ 9w -r not 



mt {At - fi cot P) Mj + |m| cos2/3sin^ 6*1^ + 



(2.59) 

where = Ml and M^ = M2. There is a difference in the sign of the off-diagonal 

tL Q tR U 

terms. This means that, as the sign conventions for the /i and A terms are the same, there 
is a difference in the relative phases of the two fields in the different conventions. Hence 
we should apply the following change to the ISAJET output: 

9t -9t, (2.60) 

i.e. change the sign of the stop mixing angle. The same argument also applies to the 
sbottom and stau mixing angles. 

We adopt the following convention for the sfermion mixing matrices: 



cos 6*^ sin 6*^ I I = ( '^^^ 

-sin (9* sin 6*^ J \qi2 J \Qri Qr2 




(2.61) 



where (jiL and qm are the left and right squark fields, for the ith quark, where i runs over 
quark flavours d, u, s, c, b, t. qn and qi2 are the squark fields for the mass eigenstates, for 
the quark i, and 6^ is the mixing angle obtained by diagonalizing the mass matrix. 

Similarly, we denote the slepton mixing as above with the matrix L^^ where i runs 
over the charged and neutral lepton flavour. /3 is the mass eigenstate and a is the left/right 
eigenstate. As we do not include the right-handed neutrino we omit the left-right mixing 
for the sneutrinos. Note that our convention for the mixing matrix is opposite to that for 
the gauginos, in the sense that here the mixing matrices act on mass eigenstates to give 
the interaction eigenstates. 

''This is the same as Ref. SPYTHIA § and SUSYGEN |]. 
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2.3.4 Higgs bosons 

The Higgs doublets are diagonalised as follows. As in Tabs. 1 and 2, Hi denote the SU{2) 
doublet interaction eigenstates which form the physical eigenstates h^,H^,A^,H^ after 
the electroweak symmetry breaking. In the charged sector: 

H+ = H+cosf3 + G+sm/3, 

= H- sin (3-G- cos /?, (2.62) 

where are the Goldstone modes. The neutral Higgs sector is parametrised in the 
interaction eigenstate as follows: 

V cos (3 + 7]i + i(i 



V2 



Hi = + (2.63) 

v2 

V is the vacuum expectation value. The diagonalisation is achieved by: 

C2 = ^°cos/3 + G°sin/3, 
Ci = ^° sin /3- GO cos/?, 
ri2 = H^ sin a + h^ cos a, 

7]i = H^ cos 0-/1° sin a. (2.64) 
Our notation is consistent with Ref. Q , but the sign of a is opposite to that of ISAJET Q . 

3. Matrix elements 
3.1 Notation 

For 2 — > 2 hard scattering subprocesses of the form 12 34, let us define the usual 
Mandelstam variables in terms of the partonic momenta as s = {pi +^2)^, t = {pi — ps)^ 
and u = {pi —pi)'^. For simplicity, the 'hat' will be omitted in the following. As initial 
state partons are taken to have zero kinematic mass, it follows that s + t + u = + m|, 
s + ts + ii4 = 0, and ut — m^m^ = spj- where = t — mg, U4 = u — m\, and px is the 
outgoing transverse momentum. 

The electromagnetic coupling is e, such that oem = e^/47r. The strong coupling is Qs, 
such that Us = 5^/47r. The QCD colour factors are Nc = 3 and Cp = {N^-l)/2Nc = 4/3. 

Colour flow directions are not stated explicitly as they are trivial in most of the pro- 



cesses whose matrix elements are given below. For a general treatment, see Ref. [13|. 
Matrix elements given here are summed over the final state colour and spin, and aver- 
aged over the initial state colour and spin. As usual, this is indicated by a long overline. 
Statistical factors of two for identical neutralinos in the final state are written explicitly. 

Let us now consider the widths T{q'^) appearing in the propagators of unstable massive 
particles with mass M{q^). In order to simplify the notation, let us define a function T{q^), 
which is defined for positive by: 

Vg^r(g2)^M(g2)f(g2). (3.1) 
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For negative q^^ the width is taken to be zero. The propagator is then given in general, 
again for positive , by: 

1 1 



(3.2) 



q2 _ M^q-i) + iy/^r{q'^) q^ - M^{q^) + iM{q-^)T{q^) 

In the HERWIG implementation, we adopt the naive running width approximation and 
neglect the running of mass, so that: 

riq"^) =q'^r{M^)/M^. (3.3) 

The following coupling notations are used in processes involving electroweak interac- 
tions: 



Qf = U{1)em charge; (3.4) 

ij = weak isospin; (3.5) 

Lf = {I^ — Qfsin^ 6w)/ cosOwsinOw', (3.6) 

i?/ = —Q/sin^ 6*^1// cos 6*^1/ sin VI/. (3.7) 

Qf and /j are as listed in Tab. 3. In terms of the above, we define the neutralino to 
sfermion coupling coefficients as: 

Lfa = {Nai cos 9w + Na2 sin 9w)Qf + {-Nai sin 9w + Na2 cos 9w) Lf, (3.8) 

Rfa = {Nal COS 9w + Na2 sin 9w) Qf + {-Nal sin 9w + Na2 COS 9w) Rf. (3.9) 

Following Ref. Q, we have neglected the contribution of the Higgsino component, which is 
proportional to the Yukawa coupling. This affects the cross sections for sfermion-Higgsino 
associated production, namely the Higgsino component of the processes gb — > tx^,bx^. 



f 


Qf 


I) 


e 


-1 


-1/2 


V 





+ 1/2 


d 


-1/3 


-1/2 


u 


+2/3 


+ 1/2 



Table 3: The C/(1)em charge Qf and the weak isospin I? for the Standard Model fermions. 



3.2 Sparton pair production 

The matrix elements and the colour decompositions for 2^2 parton-to-sparton scattering 
are given in Ref. [|l^]. The term 'sparton' refers to coloured superparticles, namely the 
squarks and the gluino. Neglecting the squark left-right mixings, the lepton-annihilation 
reaction qq* is given by: 

^ QeQq _|_ Li{Lq, Rq) 

s S-MI + iMzTz 



+ 



QiQq 



+ 



Re{Lg, Rq) 



Ml + iMzTz 



.(3.10) 
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The generalization to the case with left-right mixings is given by: 

{Lq, Rq) LgQiiQij + RqQmQRj. 
The indices i,j label the squarks that are produced. 



(3.11) 



3.3 Slepton pair production 

The matrix elements are similar to eqns. ( fj.lO ) and ( ^.11 ). For the quark annihilation 
reaction qq — > £i* one should interchange the indices {q ^ i), and change the colour factor 
from Nq to 1/Nc- In addition, for the charged current reaction qq' II^* we have 

2 

(3.12) 



\M\''{qq' ^ ^L^l 



A 2 

e sp^ 



(2sin2 Ow) Nc 



1 



s - + iMwTw 



For the lepton annihilation reaction the colour factor is 1, and there is an extra t-channel 
gaugino contribution to the production of slepton pairs of the same generation as the 
colliding lepton pair, namely. 



|A^|2(e-e^ 



eLel,eRe^) 



4 2 

e sprp 



1 



+ 



s s - Ml + iMzTz 

Re{Le, Re) 



1 

- + 



s s - Ml + iMzTz 



a=l 
a=l \ 



,0 



eLe*R,eRel) 



E 



(3.13) 



(3.14) 



The neutralino masses appearing above are the signed masses. The production of sneutrinos 
of the same generation as the colliding lepton pair is given by: 



s-M^ + iMzTz 

Re^u 



t-Ml± 



Xa 



s-M'z + iMzTz 
3.4 Gaugino production 

We consider the following seven processes in hadronic collisions: 



(3.15) 



qq 

QQ 
qq 

QQ 



Xa 'X.i 

Xi9 



(3.16) 
(3.17) 
(3.18) 
(3.19) 
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gq 



XiQ and c.c. 

xtg 

xtq and c.c. 



(3.20) 
(3.21) 
(3.22) 



c.c. refers to charge conjugation. The indices take the values a,b = 1, 2 and i,j = 1 ... 4. 
Only processes ( 3.16 ) and ( |3.17l ) are relevant to leptonic coUisions, where the colour factors 
and electroweak coupling coefficient need trivial modifications. Processes ( 3.16 ), ( 3.17] ), 
(3.18), ( 3.19| ) and ( |3.21| ) are expressed conveniently following the notation of Ref. [p!4| : 



a,f3=L,R a,f3=L,R 



(3.23) 



where, as before, Pl = (1 — 75)/2 and Pr = (1 + 75)/2. The spin-averaged amplitudes are 
given as: 



\A\^ = [\Cll\^ + \Crb.\^]u3U4 + [\Clr\^ + \CRL\^]t3U + 

+2Re[ClLCLR + C*fiiiCRL]{m3m^s). 



(3.24) 



The masses appearing above as (m^m/j^s) are the signed masses. 

The matrix elements squared for processes ( ^3.20 ) and ( [j.22 ) have the following form, 
where a labels the chirality of the final state squark: 



\M\\ 



[gsga? \ h 2(ml-mi)t3 ^ 



2Nc 



SU4 



1 + 



rrin 



+ 



M4 



(3.25) 



When left-right squark mixing is considered, the above is modified as: 

W^a ^gl — ' W?i « gh where 



2 



glQli + gRQm- 



(3.26) 



Process ( p. 16 ), qq xtXb 1 given by eqn. ( 3.241) with the prefactor {e^/Nc) and the 
following coefficients Ca/s- 



1 {SaLLq + 5aRRq){-20_ 



'/3^ 



sin 29w s- Ml + iMzT z 

i'i^ — Sai3)6quUalUbl SajS^qdValVbl 



2 sin^ 9w 

Here, following the notation of Ref. 



t-Mi 
dL 

we have defined: 



-20:1 



-K2Vb2 + 2(5a6COS Bw\ 
-Ua2Ub2 + 25 ab COS^ 6w 



(3.27) 



(3.28) 
(3.29) 



Process (|3.17|) , qq XiX'j: is given by eqn. ( [3.231) with the prefactors (e^/Nc) and 
1/(1 + 6ij), and the following coefficients Ca/3- 



Cll 



Ln 



LqiLqj 



2 ' 



sin 26*11/ s - M| + iMzT z u-M^^ 



(3.30) 
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c 



LR 



Ni3Nj3 - Nj^Nji LqjLqj 



Crl 
Crr 



sin 29yi/ 


s-Ml 


+ iMzTz 


t-Mi ' 

1L 


Rq 




- N^-iNjs 


RqiRqj 


sin 29\/i/ 


s-Ml 


+ iMzTz 


t-Mi ' 
Qr 


Rq 




- N,^Nj4 


1 RqiRqj 



sin 2ew s- Ml + iMzT z 



u ■ 



(3.31) 
(3.32) 
(3.33) 



Qr 



Process ( p. 181 ), ud — > xt)(4: is given by eqn. ( |3.23| ). The prefactors are given by 
(e'^/Nc) and (| VckmP/2 sin^ 9w) where Vckm refers to the relevant term of the Cabbibo- 
Kobayashi-Maskawa (CKM) matrix, and the coefficients Cai3 are: 



Cll 

Clr 

Crl 
Crr 



1 Ni2Val - Ni4Va2/V2 , ValL 



sm 

1 

singly 

0; 
0. 



+ 



u ■ 



s - + iMwTw 

N^2Ual+N,3Ua2/y2 

s - + iMwTw t 



Mi ' 



dL 



(3.34) 

(3.35) 

(3.36) 
(3.37) 



Process ( 3.19| ), qq — > Xi9j is given by eqn. ( |3.23| ) with the prefactor {e^glCp /Nc) and 
the following coefficients C^p: 



Cll 
Clr 
Crl 
Crr 



qi 



Mi ' 

QL 
Lqi 

^mT 

Ql 



R, 



qi 



Mi ' 
Qr 

Rqi 

-Mi ■ 

QR 



(3.38) 
(3.39) 
(3.40) 
(3.41) 



Process (|3.20| ), gq - 
Process ( ^^ , ud 



X^q, is given by eqn. (|3.25| ), with Qa = e{5aLLqi - SaRRqi)- 
> xtd^ is given by eqn. (|3.23| ). The prefactors are {e^glCp/Nc) 



and (|Vckm| /2sin 9w)i and the coefficients Cap are: 

Val 



Cll 

Clr 

Crl 
Crr 



Mi ' 

Ual 

-Mi ' 

dL 



0; 
0. 



(3.42) 

(3.43) 

(3.44) 
(3.45) 



Process ( |3.22| ), gq XIQ'l^ is given by eqn. ( |3.25D , with gL = e{5quUai + ^qdVai] 
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Processes ( 3.18 ), ( 3.21| ) and ( |3.22| ) involve flavour-changing CKM matrix contributions. 



We assume that the CKM matrix is universal and applicable also to sfermion interactions. 
Process ( ^.16 ) involves two flavour-changing vertex contributions and, in principle. 



we should sum over the internal squark flavours and external quark pairs. However, we 
believe that to a good approximation this process can be treated correctly by substituting 
an identity matrix for the CKM matrix and so this is done in HERWIG. 

3.5 Neutral Higgs production 

For hadronic collisions, the following five classes of production subprocesses are imple- 
mented: 

(a) gluon-gluon fusion via loops of heavy quarks Q = b and t and squarks Q = b and t, 

(b) ly^VF^, Z^Z^ fusion, 

(c) associated production with W^, Z^ bosons, 

(d) associated production with pairs of heavy quarks QQ, 

(e) associated production with pairs of heavy squarks QQ*. 

The scalar Higgs bosons and can be produced by all five subprocesses (a)-(e), 
whereas the pseudoscalar is only produced via the first and the last two, as the ^d^yy 
vertex, with V = W^, Z^, is absent at tree level. Higgs bosons are also produced in decays 
of heavier sparticles. 

For leptonic collisions, the following are implemented: 

(f) W^W^, Z^Z^ fusion, 

(g) associated production with Z^, 

(h) h°A° and H°A° production. 

Of these, (f) and (g) are expressed by formulae that are obtained by changing the Z^ boson 
couplings in the corresponding hadronic subprocesses. In addition, for (g), the colour factor 
1/Nc is replaced by 1. 

In the following, let us denote the neutral Higgs bosons collectively by <I>. We define 
the following notation for invariant mass squared. ^Wib] invariant mass squared 

of the particle pair [ab], Mj^^^^j = M^-^_2-^ = {pa + Pb) = s, M^afbf] = ^Pa + Pbf and 
^[a bf] ~ (P<^ ~ Pb)'^, where the subscripts i and / indicate particles in the initial and final 
states, respectively, ^^^b] negative only in the last of these three cases. 

For the gluon- fusion subprocess class (a) , the matrix elements are cast in the form seen 
in Ref. p|: 



\M\^igg^h",H") 



727rsin2 OwiN^ - l)M^ 



f 



E 



h,H Ah,H, 

q~ A~ [Tji 



(3.46) 
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f 



(3.47) 



for scalar and pseudoscalar production, respectively. The couplings gf and gj are those 
appearing in Tabs. 4 and 5, respectively. We have defined r = 4m^/M^, where M is the 
mass of the Higgs boson concerned and m is the mass of the particle inside the loop. The 
functions A are given by: 



= r + " ^^^^^^ 



■[1-t/(t)], 



for the scalars and 



(3.48) 
(3.49) 

(3.50) 



for the pseudoscalar, with the function /(r) given by 



fir) 



1 



arcsm 



log 



1 - \/r 



ITT 



T> 1, 
T < 1. 



(3.51) 



Higgs boson 


gu,i/ 




9w±,zo 


SM 


1 


1 


1 


MSSM /jO 


cos a/ sin/3 


— sin a/ cos /9 


sin(/3 — a) 




sin a/ sin /? 


cos a/ cos /3 


cos(/3 — a) 


AO 


1/ tan /3 


tan /3 






Table 4: Higgs coupling coefficients in the MSSM to isospin +1/2 and —1/2 fermions and 
to the gauge bosons W^, Z", as first appearing in eqns. ( 3.46 ) and ( 3.47] ). 



Higgs boson 


9% 


SM </> 





MSSM h° 




mjg'j + Ml cos 6 

+ (M-/(^/ni7f 


ws\n9w{LfQ\^ - i^/Q^,,)sin(a^ 
~5fa'^p)+Afmfg'j)QL^QR- 


-/?) 




f. 


'm'^fgf + A/f cos6 


>w sin OwiRfQu, ~ LfQ'i,)cos{a- 


f /3) 






















Table 5: MSSM Higgs coupling coefficients to sfermions / = g, £, as first appearing in 



eqn. ( 3.46 ). For the leptonic case, the mixing matrices Qij should be replaced by the 
corresponding slepton mixing matrices L,-,-. 
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For process (b), the matrix element for the W^W^ fusion process is given by: 

V 2 sm Ow / 



whereas for Z^Z^ fusion we have: 



I VCKM [qq"] I ' I VCKM [q'q'"] I ' M,j,j M(^,^ 

Kl3]-^W')K24]-^^) 
o2 \ 3 



(3.52) 



\MHqW2 - 939^^5) = {9%f ( . . ) (4M|) 

V4sm^ 6*14^ cos^ dw J 

The couphng coefficients gy are as given in Tab. 4. The formulae are identical to those 
for the case of leptonic collisions, with the replacement of the CKM matrix elements by 
unity for the VF^'^VF^ fusion subprocess and the substitution of the couplings L, R with the 
appropriate leptonic couplings Lg, Re for the Z^Z^ fusion subprocess. 
For process (c), the matrix elements are given by: 

(3.54) 

miqm - ^3^4) = {9tof^ ( .. „ ) ' {Ll + R-) 'fjj'^f', . 

Nc \ A sm^ Ow cos-' Ow J {s - Mlf + z 

(3.55) 

The coupling coefficients g^r arc again as given in Tab. 4. The above formula, for the Z^ 
mediated subprocess, is immediately applicable to the leptonic case with the substitution 
of the colour factor 1 to replace 1/Nc- The couplings Lq, Rq are replaced by Lg, Re in this 
case. 

For process (d), the matrix element for the qq initiated subprocess is given by: 

^2^4|'^#\2 



2sm'' OwM^s \NcJ 



(3.56) 

Here = M| for and M = Ml - 4m^ for h^,H^. The prs are the transverse 
momenta of the final state particles. These can be expressed as = 4pi -pi p2 -Pi/s — mf 
(where i = 3, 4, 5) as the incoming particles are considered massless. The propagator 
factors are defined by: 



g3 = l/(M|5]-m|), 


(3.57) 


= l/{Mf^,] - ml), 


(3.58) 


g = g3 + G4, 


(3.59) 


^5 = ^3^4/^- 


(3.60) 
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For the case Q = b, the gg initiated subprocess is evaluated with 2 ^ 1, 2 ^ 2 and 
2 — > 3 matrix elements. In order to avoid double counting, it is recommended that the 
user chooses only one of these, rather than sum over them, according to their need. For 
the 2 — > 1 subprocess the matrix element is given by: 

(3.01) 

8sm^ OwM^Nc 

The 2^2 matrix element is given by: 



2 2/ 2 2 

mi — mi I mi mi 
1 + 2^ ^ 1 + —^+ ^ 



(3.62) 



ti4 \ f3 U4 

The inconsistency in treating the initial state bottom quark as being massless and the final 
state bottom quark as being massive is unavoidable. The subprocess implementation being 
intended for high pT single 6-jet events only, the ambiguity in the treatment of the bottom 
quark mass in the hard scattering matrix element is negligible. 

We evaluated the 2^3 matrix element for gg QQ^ using helicity amplitudes as 
follows. We first wrote down the amplitude functions using the corresponding Feynman 



rules. We then adopted the convention of Ref. |16| for the heavy quark helicity eigenstates, 
and adopted the circular polarization for the gluons. This allowed some simplification in 
the computation which, together with the factorization of amplitudes, made the code very 
much faster than any other code which is publically available. 

For process (e), the matrix element for the qq initated subprocess is given by: 



\M\^{qiq2^Q3Ql^t 




sin^ 9wM^s 

X [IGsWts + IG^Wta - 2Re {giGi)pT3 ■ PTi] ■ (3.63) 
The propagator factors are given by: 

G3 = 1/ {{P4+P5? -ml + imsFs) , (3.64) 
= 1/ {{ps + Pbf -m\ + imiVi) . (3.65) 
The dot product of two transverse momenta is two dimensional, i.e. pTs ■ Pta = PxsPx^ + 
Py^Py-i- 

The gg initiated subprocess is again calculated using helicity amplitudes in an opti- 
mized basis. This time the two orthogonal polarisation vectors for the initial state gluons 
are chosen to be in the x and y directions, yielding a relatively compact expression: 



\M\Hg1g2 ^ Q3Ql^5) 



4sin2 9wM^{Nl; - 1) 

Y. E 

(3.66) 



i=x,y j=x,y 



\At{i,j)\^ + |A(i,j)P - -^\At{i,j)+Au{i,j] 
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The colour flow components are distributed using the method of Ref. |13], namely, the 
{1/Nq) term is distributed amongst the t-channel and u-channel colour flow components 
by their ratio. The two amplitude functions At and Au are given by: 

= 2 (pi3Pj4^24^13 +Pi3Pi3^13^45 +Pi4Pi4^24^35) 

- ^ (^23^45 - Pz^Gsb) - ^ (^35 + ^45) , (3.67) 
•^u{i,j) = 2 (pi4Pj3^l4fe +Pi3Pj3fe^45 +Pj4Pi4^14fe) 

+ ^ (^,3^45 - PZ4Q35) - {G35 + Ga5) . (3.68) 
\/S z 



The propagator factors are given by: 



G13 


= 1/{{P1+P3? 


-ml 




G23 


= 1/((P2+P3)' 


-ml 






= l/{{Pl+PAf 


2 

- m^ 




Q2A 


= ^/{{P2+PA? 


- ml 




^35 


= 1/((P3+P5)^ 


- m\ 


+ zm4r4) 


^45 


= 1/((P4+P5)' 


-ml 


+ im3r3) 



(3.69) 
(3.70) 
(3.71) 
(3.72) 
(3.73) 
(3.74) 



For process (h), the matrix element is given by: 



Ll {Ll + Rl) (cos2(/3 - a), sin2(/3 - a)) 



{s - Mlf + MlTz 

We have used the same notation as in eqn. ( p.lOp . The coupling factor cos^(/3 — a) is for 
h^A^ production whereas sin^(/3 — a) is for H^A^ production. 

3.6 Charged Higgs production 

For hadronic collisions, the following production modes are implemented: 

(3.76) 
(3.77) 
(3.78) 
(3.79) 
(3.80) 
(3.81) 
(3.82) 
(3.83) 

Charged Higgs bosons are also produced through the decay of heavier objects including 
gauginos and third generation sfermions. Perhaps most notably, an important potential 
production mode for charged Higgs bosons is through the decay of top quarks and, as 



gh- 


-> tH^ + c.c; 


99 - 


->■ tbH~ + c.c; 


qq - 


->■ tbH~ + c.c; 


99 - 


^tb*H' + c.c. 


qq - 


■^tb*H- + c.c. 


qq - 




bb - 




qb- 


-> q'bH'^ + c.c. 
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with other decay processes, this is achieved by specifying this decay mode in the input 
file. The decay is carried out using the three-body matrix element. This matrix element 
is not included in the spin correlation algorithm |]l^ which is the default treatment from 
version 6.4. If this decay is required, the matrix element code NME=200 should be used 
and spin correlations and super symmetric three body decays should be switched off, i.e. 
SYSPIN= . FALSE . and THREEB= . FALSE . . 

When simulating processes ( 3.77 ), ( 3.78D and ( p. 83 ) one should bear in mind that there 
is potentially a large contribution from top quark decay when this decay mode is available. 
In this case the user should choose to simulate either the top quark production or the 
relevant three-body process in order to avoid double counting. In the limit in which the 
bottom quark is generated by the collinear splitting g — > bb, process ( 3.77] ) is described 
by the two body process ( 3.76| ). This contribution becomes significant if the mass of the 
charged Higgs boson is comparable or greater than the top quark mass. Again, the user 
should choose to simulate either of the two processes in order to avoid double counting. 

For leptonic collisions, in addition to the top quark decay, charged Higgs bosons can 
be produced in a process analogous to process ( p. 81 ): 



e e 



(3.84) 



The matrix element is obtained immediately from the hadronic case by changing the ap- 
propriate photon and couplings and removing the colour factor I/Nq. 



Process (3.76) is given by: 



\2Nc) 



|VbKM[t6]P 



2 f ml tan^ P + cot^ (3 



2M^ 



-Ua 



1 + 2- 



li4 



rriA 



1 + -^ + — i 



is 



li4 



(3.85) 



Process ( |3.77 ) is evaluated using helicity amplitudes in an optimized basis, similarly 



to the evaluation of the neutral Higgs boson production process gg —>■ QQ^ in Sect. 3.5 



Process ( p. 78 ) is evaluated using a more general formula for qiq2 QsQi^ which is 
given by: 



+ 2 (p3 .p4-Am3m4) [-s\g\^ + 2Re {gg*spl, + GGIp^Ts " GsGlPn)] 



-4Re 



GG3PT4 + GGIPt. + ^ {S3P0 • Pa + QIpo • P3, 



pk+Ph-PTs 



Aim [GGlPzi - QQlPzi] A'\/s ijPxiPyj, - Px-iPvi 



(3.86) 



with pq = Pi + P2- Denoting the QsQ^^ vertex as {gs + 9P75) 
coefficients appearing above are defined by: 



the various coupling 



9p, 



(3.87) 
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~ 12 ■ 



X' 



(3.88) 
(3.89) 



Process ( |3.79| ) is evaluated using helicity amplitudes. The expressions are the same as 
eqns. ( |3.66D -( ^.74D but the coupling coefficient needs to be changed. Let us consider the 
vertex involving an incoming and outgoing tibj. We have the following replacement: 



1 

71 



sin 2/3 



m? tan /? + cot /? 



+ Ab tan P) Q\,Q% - ^ (-^ + At cot /?) Q'^.Q^^. 



(tan /3 + cot /3) Q^Qni 



wife 



(3.90) 



Process ( 3.81 ) is given by: 



A 1 



Nc 

QqQH+ 



+ 



s- + iMzTz 



+ 



QqQ 



+ 



RqLH+ 



M| + iMzTz 



.(3.91) 



We have used the same notation as in eqn. ( p. 10 ). Qh+ = +1) and Lh+ is given by 



eqn. (|3.6D with /| = +1/2. For the lept onic case the subscripts q should be replaced with 
e, and the colour factor 1/Nc should be replaced by 1. 

The matrix element for process ( p. 82 ) is as given in [p^ ]. 



Process ( p. 83]) is given in [19|. We use a more compact expression which sets the 
kinematic bottom quark mass to zero. For the process hiU2 h^D^^H'^ we have: 



\M\''{hiU2^hD^H-^ 



Re(2^4^t 



1 



{Ml 



[24] 



\VcKM[th]\'W + l^-f P + Re(VcKM[t6] • 2AlAt, 

2 



/mfot^y [3545] ^ 



mw I 



,2 , ^2^, ,.2 ^[2545]. 



(e',H + ei)(-^[i3] 

""'^^ (^h.H + ^a) \2Re{Qt){Ml^. - Ml 



(3.92) 
(3.93) 

(3.94) 



-Re(gt)M25][1243] +Im(gt)M55j •4det(1234) 



(3.95) 



Here we have defined [abed] = 4{pa -pb Pc-Pd+Pa -Pd Pb -Pc-Pa-PcPb-Pd) and det (a6cd) 



Hay pa ^ iJiV paPaPlPcP^ whcrC €0123 = +1- 
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The propagators are defined by: 



Qt 

Qh,H 

Qa 



^[35] 



imtVt 



TUf) sin a cos(/3 — a)/Mw cos /? 



2\2 



^2 ' 



+ 



TJib COS a sm{(3 — a)/Mw cos /? 



^[13] 



tan /3/Mi 



^[13] 



Ml 



Processes with antiquarks in the initial state can be obtained by exchanging the subscripts 
(1 <-> 3) for the subprocess bU — > bDH^, (2 <^ 4) for the subprocess bD bUH~^, and 
both for the subprocess bD bUH^. Note that the definitions of the invariant masses 
squared are such that these imply momentum substitutions of the form (pi <^ —ps)- 

3.7 R-parity violating processes 

The cross sections for the implemented R-parity violating hadron-hadron processes are 
presented in Ref. [20|. In hadron-hadron collisions only those 2^2 processes that can 
occur via a resonant s-channel sparticle exchange are included, i.e. processes that can only 
occur via a t-channel diagram or for which the s-channel resonance is not kinematically 
accessible are neglected. 

In e^e~ collisions we include all the possible 2^2 single sparticle production mecha- 
nisms, even if there is no s-channel resonance, and a limited set of 2 — > 2 processes where 
the final state contains two SM particles. As with the hadron-hadron processes, there 
are four types of process: single gaugino production; single slepton production in associa- 
tion with a gauge boson; single slepton production in association with a Higgs boson; SM 
particle production. 

3.7.1 Single gaugino production 

In e~^e~ collisions, there are two possible production mechanisms for a gaugino and a 
lepton, either e'^e~ or e^e~ X^^~ ■ The matrix elements for these processes 

have been calculated previously in Refs. [^, |2^. We have recalculated these matrix 
elements and the results are in agreement with 



We also agree with Refs. |2l|, p2| , 
which do not include sfermion mixing, if we set our mixing to zero. All matrix elements 
are averaged over the spins of the incoming leptons as before. 
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+ 



i{v*)a{li2)R{yj) [s- m-p^ ] st 



(3.96) 



R{u,)s\{\a{u*)\' + \b{u*r) 



9 \ji 



s-mj -ML 

'^J Xi 



J Xi 



+- 



, , ML 

2 V X, 



u) (ml-u) [\am' + m)\'] 



2s 



U-M3 



Mi 



(3.97) 



The couplings a, h of the sleptons to the gauginos are given in Ref. ||2^, and 

R{a) = \ 5 . 

{s-Mlf+TlMl 

3.7.2 Slepton production in association with a gauge boson 



(3.98) 



There are a number of production processes for a slepton in association with a gauge 
boson. We have confirmed the cross section for these processes, which were first calcu- 
lated in Ref. ||2^. The sneutrino resonance is not usually accessible for these processes 
as the charged slepton-sneutrino mass difference is less than the mass. However, the 
resonance can be accessible for the third generation due to left -right mixing of the staus. 

The matrix elements for these processes are given below, including left -right sfermion 
mixing: 
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Ml 



Ml 



2M| 



Ml 



(3.101) 



where pcm is the momentum of the final-state particles in the centre~of-mass frame. As we 
neglect the electron mass, the cross section for e~e^ — > ^Vj is divergent as u, t — > 0, i.e. in 
the limit that the photon is collinear with the incoming lepton beams, and therefore a cut 
on the transverse momenta of the outgoing particles is used to regularize the divergence 
for this process. The couplings of the to the leptons and sleptons are given in Ref. pO|. 



3.7.3 Slepton production in association with a Higgs boson 

There are a number of processes in which a Higgs boson can be produced in association with 
a slepton. The cross sections for these processes tend to be small due to the small lepton 
masses and the inaccessibility of the sneutrino resonance as the charged slepton-sneutrino 
mass difference is smaller than the charged Higgs mass. 
The matrix elements are given below 



|A^|2(e: 



2j-l|2 



\Mnel 



V*Hf) 



\E: 



l\2 



ut 



ut 



+ 



\E[ 



i\2 



ut-MlM^o 



(3.102) 



. (3.103) 



The couplings H of the Higgs bosons to the sleptons are given in Ref. ||2^ and the couplings 
E\ of the leptons to the Higgs bosons can be obtained from the couplings given in Ref. |20| 
for the down- type quarks by replacing the quark mass with the lepton mass. 



3.7.4 SM particle production by sparticle exchange 

As with the hadron-hadron processes described in Ref. i2' 



we have only included 2 — > 
2 processes where two SM particles can be produced by a resonant s-channel exchange. 
Where this is possible we have included any f-channel sparticle exchanges, and the relevant 
interference terms. We have also included the relevant SM diagrams and the interference 
between the SM and the R-parity violating diagrams. We have generalized the results of 



Ref. ||2J, giving 



^ {(1 + co,ef [|/£^|2 + + + ^ 2Re(/r^/L) + 2Re(/^l/k)] 

+ (1 - COS0)2 + + 4 + |^|^^|2] \ (3_104^ 
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where the coefficients are given by 
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j=i,3 2e2 [s-Ml^^iV^M^^ 
for e~ el — > dfcd/. The colour factor Cfac is 1 for charged lepton production and Nc for 



quark production. The left/right charges of the fermions are given in Sect. 3.1. This 
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formula corrects the sign of the interference term given in [p4|| , which is necessary both 
to reproduce their numerical results and to obtain the correct SM Bhabha cross section 
formula in the limit that the R-parity violating couplings are zero. 



3.8 Decay matrix elements 

The decay matrix elements are implemented for both R-parity conserving and R-parity 
violating three-body decays of the gluino and the gauginos. The R-parity conserving decays 
are performed using the helicity amplitude expressions given in Ref. [17|. The R-parity 
violating formulae are as given in Ref. pO| . 

Spin correlations in the distribution of cascade decay products are implemented in all 
decay processes, as well as the R-parity conserving 2 — > 2 SUSY production processes. The 
implementation is as described in Ref. |17]. 



4. Conclusions 

The Monte Carlo event generator HERWIG has been extended to include the simulation of 
supersymmetry. We have summarised its salient features. We have defined the parameter 
conventions and listed the matrix elements for all original calculations. 
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A. Process codes 



We summarize the process codes IPROC for MSSM hard scattering processes implemented 
in HERWIG. 



IPROC 



700-99 
700 
710 

706+4IN1+IN2 

730 

728+2IC1+IC2 

740 
736+5IL 
737+5IL 
738+5IL 
739-t-5IL 
740+5IL 

760 
757+410 
758+4IQ 
759+4IQ 
760+4iq 



800-99 

800 

810 
810-hIN 

820 
820+IC 

830 

840 

850 

860 

870 

867-h3lLl-ML2 
880 

877+3IQ1+IQ2 



910 
920 



960 
970 



955 
965 
965 



Process 



R-parity conserving SUSY processes 

^ 2-sparticle processes (sum of 710-760) 

neutralino pairs (all neutralinos) 
£'^£~ XiNiXiN2 (INI ,2=neutralino mass eigenstate) 
£~^£~ chargino pairs (all charginos) 
£~^£~ XiciXic2 (lCl,2=chargino mass eigenstate) 
£'^£~ slepton pairs (all flavours) 
£+£- ^ £l,r£%r (IL = 1, 2, 3 for £ = e, jl, f) 
£+£- £l£*l (IL as above) 
£+£- £l£*r (IL as above) 
£+£- IrI*i^ (IL as above) 
£+£- VlvI (IL = 1, 2, 3 for V^, V^) 
t^l~ squark pairs (all flavours) 
£+£- ^ QL.Rqi.R (IQ = 1-6 for q = d.J) 
£~^£~ qLql (IQ as above) 
£'^£~ Ql^h (IQ as above) 

qRq*R (IQ as above) 



R-parity violating SUSY processes 

Single sparticle production, sum of 810-840 



£+t 
£+t 
£+£- 
£+£- 

e+t 

£+t 
£+£- 



(all neutralinos) 
Xjj^z^i, (IN=neutralino mass state) 
XTel (all charginos) 
Xic^t ^ (IC=chargino mass state) 

v.Z"^ and l+£- F+VK" 
D,hPlH^IA° and £+£- IfH- 

Sum of 870 and 880 

e+e- ^ £+£-, via LLE only 

£+£- ^J^i^fL2 (ILl, 2=1,2,3 for e,n,T) 

£+£- dd, via LLE and LQD 

£+£- diLidiL2 (IQl, 2=1,2,3 for d,s,b) 



£+£- 
£+£- 



£+£- 
£+£- 



£+£- 
£+£- 
£+£- 



H+H- 



Table 6: The MSSM hard scattering processes implemented in HERWIG (continued on following 

pages). 
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IPROC 



Process 



R-parity conserving SUSY processes 

2-parton 2-sparticle processes (sum of those below) 

2-parton 2-sparton processes 

2-parton — > 2-gaugino processes 

2-parton — > 2-slepton processes 



3000-999 
3000 
3010 
3020 
3030 



3100+ISQ 



gg/m-^ m {cf. Rcf. |1| for definition of ISQ) 



3200+ISQ 



gg/qq~>qq*h,H,A (") 



3310,3315 
3320,3325 

3335 

3350 

3355 
3360,3365 
3370,3375 



qq 
qq' 
qq' 
qq 
qq 
qq 
qq 



Z"H°A°H" 



ft," + ch. conj. 
b + ch. conj. 

A" + ch. conj. 
t + ch. conj. 



3410 
3420 
3430 
3450 



hg- 

bg 

bg- 

bg 



3500 



bq bq'H^ + ch. conj . 



3610 
3620 
3630 



qq/gg 

qq/gg h" 
qq/gg -4° 



3710 
3720 



qq 
qq 



q'q'h° 
q'q'H" 



3810+IQ gg + qq ^ QQh° (IQ for Q flavour) 

3820-MQ gg + qq^QQH° {") 

3830+IQ gg + qq~^QQA"C) 

3839 gg + qq btH^ + ch. conjg. 

3840-l-IQ gg -> QQh° (IQ as above) 

38504-IQ gg^QQH° (") 

3860-1-iq gg->QQA^ (") 

3869 gg btH+ + ch. conjg. 

3870-Mq qq QQhP (IQ as above) 

3880-MQ qq-^QQH^C) 

3890+IQ gg->Qg^°(") 

3899 qq -> btH+ + ch. conjg. 



Table 6 continued. 
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IPROC 


Process 


4000-99 


R-parity violating SUSY processes via LQD 


4000 


single sparticle production, sum of 4010-4050 


4010 


Ujdk 


x'^li , djdk (all neutralinos) 


4010+IN 


Ujdk 


Xjj^Z", djdk XiN'^i (lN=neutralino mass state) 


4020 


Ujdk 


djdk X~(^t charginos) 


4020+IC 


Ujdk 


Xic>^i^ djdk — ^ Xic^t (lC=chargino mass state) 


4040 


Ujdk — » 


f+ZO, Ujdk 'i>^W+ and djdk ^ 


4050 


Ujdk 


I+h^/H^/A"^, Ujdk ViH+ and djdk I+H- 


4060 


Sum of 4070 and 4080 


4070 


Ujdk — > 


uidm and djdk ~* didm, via LQD only 


4080 


Ujdk 


I'jl'i^ and djdk ~> ^t^k ^ '^i^ LQD and LLE 


4100-99 


R-parity violating SUSY processes via UDD 


4100 


single sparticle production, sum of 4110-4150 


4110 


Uidj 


x'^dk, djdk X^Ui (all neutralinos) 


4110 -fIN 


Uidj 


XiN^fc, djdk XmUi(lN as above) 


4120 


Uidj 


X^Uk, djdk x~di (all charginos) 


4120 +IC 


Uidj 


Xj^Ufe, djdk Xicdi (IC as above) 


4130 


Uidj 


gdk, djdk -> gui 


4140 


Uidj — !■ 


hlZ°, djdk ilZ°, Uidj i*W+ and djdk -> b*W- 


4150 


Uidj — > 


dl^h'^/H^/A^ djdk ^ u*^hyHyA°, u,dj ^ ul^H+, djdk ^ d*^H- 


4160 


Uidj — > 


uidm, djdk — > didrn via UDD. 



Table 6 concluded. 
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